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Abstract 

QJ , We propose and investigate the thermodynamic Bethe ansatz equations for the min- 



imal models (associated with the A^^i Lie algebra) perturbed by the least (Z 



N 



^ , invariant) primary field <I>jv. Our results reproduce the expected ultraviolet and in- 
' frared regimes. In particular for the positive sign of the perturbation our equations 



describe the behaviour of the ground state flowing from the model to the next 

P- 



W^_i fixed point. 
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Recently Al.Zamolodchikov has pointed out the importance of the thermodynamic Bethe 
ansatz (TBA) in the context of the non-diagonal scattering theories [1,2]. In particular, 
for the RSOS models, the respective TBA equations for the ground state reproduce the 
expected ultraviolet and infrared behaviours. The RSOS scattering theory describe [1,3] 
the scaling behaviour of the minimal conformal Mp field theories perturbed by the 0i 3 
operator. For the negative sign of the perturbation, the TBA equations describe the massive 
behaviour of the (p-1) Z2 symmetric kinks. On the other hand, for the positive sign, they 
describe the Renormalization Group (RG) trajectory flowing from the Mp critical point 
to the Mp_i fixed point. This is a well known example of RG flow between two non- 
trivial fixed points [4,5]. It seems important to generalize this approach to other systems 
possessing "higher" symmetries. The natural candidates are the Z^r-invariant conformal field 
theories related to the minimal models [6] with central charge c— [N — — ^p^^i^ ) , 
p = A'"-|-l,A^-|-2,... . The analogue of the (/>i,3 operator in the theories is the least (Zjv 
invariant) relevant primary field associated with the weight of the adjoint representation 
of SU{N)p^i_N, which has conformal dimension = 1 — For the sake of simphcity 
we first consider the Wp theory. The perturbed action is defined by 
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Aw3 = Ac FT j 



where the field $3 is identified with the operator 



1 2 

1 2 
1 2 



in the notation of ref.[6]. 



The field 



1 2 

1 2 



preserves the global symmetry of the Wp model and therefore the 



perturbed action (1) is still Z3 invariant. Moreover, the action (1) defines an intcgrable 

theory [7,8,9] and its physical behaviour depends on the sign of the coupling constant A. 

For A < the action (1) describes the scaling region of a generalized Ajv_i-RSOS [10] model, 

and its kink-kink S-matrices have been proposed in ref.[ll] . In the case of A > 0, based 

1 2 



on the product expansion properties of the field 



1 2 



and the RG analysis [8,12], 
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one expects that the theory flows to the Wp_i fixed point [8,12,13], p > 4. We will 
henceforth refer to these theories as VF^ '■^^ models. Up to now there has not been any general 
scheme in order to derive the TBA equations directly from the non-diagonal S-matrices. 
However, Al. Zamolodchikov has pointed out that (see also refs.[20],[21]) in the usual RSOS- 
models [1,2,19] , the final result assumes a general form in terms of coupled integrals equations 
with a universal kernel, respecting the symmetries of the kink interactions. Then, using the 
available symmetries and the knowledge of certain known limits it is not difficult to guess 
the associated TBA equations. The vahdity of these equations can be confirmed through 
checks done a posteriori. For example, analyzing the ultraviolet and infrared behaviour of 
the proposed TBA equations and also comparing them with the respective perturbation 
theory. The TBA equations that we propose for the action (1) are the following 

em - KiO) + ^ E kj E / d^'KaAO - 0')Ln{l + e-^'^''^) = (2) 

where kj = 5ij - 5ij+i - 6ij^i and Kij{x) = ^2,2(2:) = - 2cosh{x)-i '^ ^1,2(2:) = K2,i{x) = 
-2^Jkri- = ^imRcosh{9) (A < 0) and i.^ = ^S,,, + ^5,,,_3 (A > 0). 
The corresponding free-energy E{R) at temperature T — ^ is given by 

= / deu^M^ + e-^'^^'^) (3) 

i=i b=i-' 

The Zs-symmetry is expressed by the kernel Ka,}, properties and by the fact that e\ — e^, 
i — l,2,...,p — 3, is a solution of Eq.(2). For p=4 we automatically recover the TBA 
equations of the 3-state Potts model perturbed by the thermal field 02,i = § [14] . 

The first check consists of an analysis of Eqs.(2,3) in the ultraviolet limit, R ^ 0. In 
this limit the frcc-cnergy E{R) is expected to behave as E{R) ~ — ^ [15], where c is 
the associated central charge. Using the standard calculations of the leading behaviour 
i? — > of the free-energy E(R) [14,16], the central charge can be written in terms of the 
Rogers- Dilogarithm function L(x) as 



p-3 - 

cuv = 2 E 

i=l . 



(4) 



and 

^ /-a ^ Ln{x) ^ Ln{l — x) 
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L(a) = ^ dx 



1 — X X 

where the constants Xi — e~^<^°^ = e"^?*^^^ and yi = e~'^i^°°^ — 6"^**^°°^ are given by 



^« + l = n ■ /t\ a=l,2,...,p-3 (6) 

and 

yi = 0; y„+l=n wf^' ^ = 2, 3, 3 (7) 

i=i p j 

Finally, using the Rogers Dilogarithm sum rule 

\^T(^^\ (p-2)(p-3) 

we find the central charge cuv of the minimal Wp models 

=- = 2(l-^)> f = i5,... (9) 

Turning to the infrared limit i? — > oo, for the Wp^~\ it is straightforward to recover 
the massive behaviour c — 0. In the case of Wp^'^\ however, one may notice that the 
pseudoenergies ei (e^-s) behave as {^^e~^) in the region —Ln{MR) < 9 < Ln{MR). 

Proceeding as in the ultraviolet regime, our final result is 

= 2 - -^ii^j , p = 5,6,... (10) 

as it is expected for the RG fiow 

Now we start to discuss the next to leading order in R corrections to the free-energy E(R). 
Using general properties of Eq.(2) [17], it is possible to show that the functions Y-"-{9) = 
e~^?(^) have a Laurent expansion in terms of the variable t — ep+^ . As a consequence [17] we 
expect that the scahng function F{r) — ^^f^ (r=MR or r=mR) has a series expansion in 

6 „ 

power of , besides a usual bulk term proportional to R . Moreover, when p = 3i — 1, i — 
2,3, ... we also expect a logarithm divergence at order ^±1^ which gives a term of r^Ln(r) 
form. 
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For the sake of clarity, here we discuss in details the simplest case of models. 
From our considerations above, we have the following expansion for F(r) 



Fir) = -^-^,r'Ln{r) + Y.frr\ W^^ (11) 



10 Stt 



i=2 

F(r) = -—- - -1-r^Lnir) + V /+r\ 14/5+^ (12) 

^ ' 10 247r 87r2 ^ ' 'H,'^' ' ^ ^ ' 

The bulk term in the W^^^^ comes from a simple manipulation of Eqs.(2,3) [2,14]. 
However, the calculation of the constant A of the j^^r"^ Ln{r) term is a bit more involved. 
Adapting the same method used by Al. Zamolodchikov for the RSOS model [1,2], we find 
that A = — I {! +—) = — "ttQ- The coefficients and /" can be calculated solving 

numerically the Eqs.(2,3). In table (1) we collected the first 5 of them. We notice (within 
our numerical precision) that f2 = f2 ) then we have the same mass scales (m=M) for 
both signs of the perturbation A. It is also important to compare Eqs. (11) and (12) with 
the expected results from the perturbation theory for the VF5 models. Comparing the 
logarithm divergence in the perturbation theory with the term —-^r'^Ln{r), we find the 
exactly relation between the mass m (M) and the coupling constant A 

m = M = (27r)A (13) 

and from the standard results of perturbation theory [1], the coefficients ff are given by; 

ft = ^(^)'^*3,*3,*37(^)' (14) 

where C$3,#3,#3 = |7( i^)7(^)7(f )^0, and7(x) = We estimate fi' = t3.71575994£;- 

03, in accordance with the TEA results from Table 1. Moreover, we observe that the in- 
terchange of signs of form /+ = (— 1)"/~, ??, = 2,3, ... (see Table 1) is in agreement with 
what one expects from the perturbation theory around the positive and negative signs of A. 

^For W^p''^'' theory, the general form of the r^Ln{r) term is (— 1)^" -nrA^r^Ln(r),p = 5,8, ... 

^This structure constant may be calculated using the decomposition = | = 5 + ^ [22] 

•^We have also checked that, for the next model Wq, the coefficients ff^ are in agreement with the 

perturbation theory. The corresponding structure constant may be calculated noticing that c-^^^ — ^ 

appears in the superconformal minimal models SMp, p = 12. [23] 
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In order to get some insight about the next to leading corrections in the infrared regime, 
we have solved numerically the Eqs.(2,3) for large R, namely in the region 50 < i? < 400. 
The first correction is determined to be with very high precision. This suggests that 
the infrared behaviour of the W^^~^^ model may be seen as the theory perturbed by the 
field TT (plus other irrelevant fields) 

A 3(+) = Ay^-i + / (fxTT{x) (15) 



where T{T) is the holomorphic (anti-holomorphic) component of the stress-energy tensor of 
the Wi model. 

We stress that a similar situation happens in the tricritical Ising model flowing to the 
Ising theory [1,2]. We can check this last proposal as follows. From perturbation theory [1,2] 

2 

the coefficients ai,a2 proportional to the terms aiR~^ and a2-R~^ satisfy the relation ^ = 

2 

— ^ = — ^. Our numerical calculation predicts ^ = — 0.033332(±1), in good accordance 
with the perturbation theory. 

For general jo > 5, the infrared behaviour is mainly governed (as discussed before) by 
the pseudoenergies e] = ef,i = 2, 3, ...,p — 4. This changes the periodicity of the function 
Yl^{9) = e~'^i^^\ and predicts a new expansion to F(r) in powers of R~~. This exponent 

2 1 



agrees with the conformal dimension of the 



operator in the VFi_^ minimal models. 



2 1 

Then, we may argue that this field is one of those that are responsible^ for the infrared 
corrections in the W^^^^ models. 

The generalization of our results for arbitrary N is as follows. The TBA equations are 
written as 

1 p-N N-l 

em - KiO) + ^ E h E / d(^'Kaii0 - e')Ln{l + e-^?(^')) = (16) 
and the Kernel K^'j, is given by 

KZ{e) = -i^^LniSaM): Km = ^-[^Ln{Sa,m, ^ + J (17) 



*Such argument has been originally used by Al.Zamolodchikov in the case of the usual RSOS models [1,2] 



where Sa,b is the two-body S- matrix of the Z(N)-parafermionic field theory perturbed by 
the energy operator with anomalous dimension A — [18]. The Sa,b has the following 
expression 

rnin(a,b) Sh-{6 iZLHl^ 

Sa,b{0) = n F\a-b\+2l-l{0), F{e)m = QA(n,i^m\ (l^) 

One can notice that the only poles in the physical strip are due to the two-body S- 
matrices Sa,boi the Z(N)-model. Then, the functions have the form 

= 5i,rmaRcosh{e), (A < 0); < = ^e^5,,i + ^e-^5,,,_Ar, (A > 0) (19) 

where the masses rua and Ma are given in terms of the masses ratios of the Z(N) theory, 
namely m„ = rn'-^\ ; M„ = M^^, a = 1, 2, N -I. The free-energy E(R) keeps the 
same form as Eq.(3), only substituting Y^^Zi by the general form Y^iZi ■ 

As in A'" = 3 case, the Z(N)-symmetry of the pseudoenegies (j = 1,2,..., A — 1) 
are guaranteed by the properties of the kernel K^^^- We also have two particular limits of 
Eqs.(16,17) that are easy to be taken. For N — 2, the Eqs. (16,17) are the same ones proposed 
by Al.Zamolodchikov to describe the minimal models perturbed by the 01^3 field. In the 
case of arbitrary N and p = A^ -|- 1, they reproduce the known results for the Z(N)-models 
perturbed by the leading energy operator [14,18]. 

Proceeding as in the case A^ = 3, we have checked the ultraviolet and infrared behaviours 
of Eqs. (16, 17). We find that both regimes are in accordance with the central charge of the 
minimal models (p=N-|-l,N-|-2,...) related to the VFAjv-i- algebras. 

In summary we have proposed and investigated the TBA equations of the minimal mod- 
els perturbed by the least {Z^ invariant) primary relevant field The introduction 
of an extra degree of freedom, i.e. the "fiavour" a = 1, 2, A— 1, is the basic characteristic 
of Eqs. [2, 16], which allows one to generalize our results for other symmetries. In our con- 
siderations we have presented only the final results without entering into technical details. 
Following ref.[19] it is also possible to extend our results to a general class of coset models 
based on the SU(N) group. A detailed account of the technical part of this letter as well as 
the generalization to the SU(N) coset will be published elsewhere [24]. 
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Tables 

Table 1- The coefficients ff, and for 1=2,3,4,5,6 



i 


/r 




2 


2.92560874(±l)^- 02 


2.92560875(±l)^-02 


3 


-3.7157599(±1)E - 03 


3.7157603(±2)£; - 03 


1 


1.68716(±2)E - 01 


1.68712(±2)E - 01 


5 


1.4368(±1)£; - 04 


-1.4367(±l)E-04 


6 


1.645(±1)E- 05 


1.642(±1)£; - 05 



9 



